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Abstract. Although controversial, the scenario of microlensing as the dominant mechanism for the long-term optical variability 
of quasars does provide a natural explanation for both the statistical symmetry, achromaticity and lack of cosmological time 
dilation in quasar light curves. Here, we investigate to what extent dark matter populations of compact objects allowed in 
the currently favored Sl M = 0.3, C1 A = 0.7 cosmology really can explain the quantitative statistical features of the observed 
variability. We find that microlensing reasonably well reproduces the average structure function of quasars, but fails to explain 
both the high fraction of objects with amplitudes higher than 0.35 magnitudes and the mean amplitudes observed at redshifts 
below one. Even though microlensing may still contribute to the long-term optical variability at some level, another significant 
mechanism must also be involved. This severely complicates the task of using light-curve statistics from quasars which are not 
multiply imaged to isolate properties of any cosmologically significant population of compact objects which may in fact be 
present. 
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1. Introduction 

In principle, the possibility that the long-term optical variability 
of quasars may be due to microlensing caused by a population 
of planetary-mass compact objects (Hawkins I1993I Hawkins 
I1996I has several attractive features. It provides a natural ex- 
planation for the symmetrical and largely achromatic nature of 
quasar variability on time scales of a few years, in addition to 
explaining the lack of cosmological time dilation in the vari- 
ability time scales (Hawkins 2001 1. This idea is nonetheless 
controversial, since it appears to require a cosmologically sig- 
nificant population of dark, compact bodies in a mass range 
which due to constraints from the MACHO and EROS projects 
(e.g. Alcock et al. 1998 Lasserre et al. 2000 1 can contribute 
no more than a small mass fraction to the dark halo of the 
Milky Way. Inferring a limit on the cosmological density of 
compact objects from the mass fraction of a single halo is, how- 
ever, difficult given our poor knowledge of the cosmic variance 
of dark matter properties. There are also concerns that more 
diffuse lenses (e.g. gas clouds) would not produce sufficiently 
strong microlensing effects to be detectable in our own halo, 
but only at cosmological distances (Walker & Wardle 1999 1. 
Gibson & Schild ( 1999 1 have furthermore suggested that strong 
clustering of lenses inside our halo could make a large popula- 
tion of planetary-mass bodies evade the MACHO/EROS con- 
straints, yet allow detection of such objects in studies of quasar 
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microlensing, where the optical depth due to lensing is much 
higher. 

Several papers critical of the scenario of microlensing- 
induced quasar variability have been published. Lacey {l994i 
& Alexander ( 1995 1 both modelled the relation between time 
scale and redshift predicted by microlensing scenarios and 
found the result to be in conflict with observations of Hawkins 
( 1993). The observed relation has since been slightly revised 
(Hawkins 1996 ), and neither Alexander nor Lacey explored the 
entire parameter space relevant for the microlensing scenario 
or the selection effects involved. Baganoff & Malkan (1995 1 
argued that the observed redshift trend could be explained by 
intrinsic, wavelength-dependent variability, but this was later 
refuted by Hawkins & Taylor ( 1997 ) on the basis of predicting 
color-dependent variation time scales in conflict with observa- 
tions. 

Despite the qualitative agreement between statistical prop- 
erties of long-term optical quasar variability and expectations 
from microlensing scenarios, the detailed comparisons per- 
formed between observations and microlensing simulations 
have so far been very limited. Here, we use numerical simu- 
lations of quasar microlensing to investigate the full parameter 
space of relevance for the currently favored Qm = 0.3, Qa = 
0.7 cosmology and investigate to what extent microlensing re- 
ally can reproduce quantitative features of the observed-long 
term variability, like the average structure function, the ampli- 
tude distribution and the amplitude-redshift relation. 
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In Sect. 2, the machinery used to simulate the microlensing- 
induced optical variability of quasars on long time scales is out- 
lined. Section 3 describes the parameter space explored and the 
method of building synthetic quasar samples used in compar- 
ison with observations. In Sect. 4, 5 and 6, the observed aver- 
age structure function, cumulative amplitude distribution and 
amplitude-redshift relation are compared to the corresponding 
results from our microlensing simulations. In Sect. 7 we dis- 
cuss various assumptions on which this investigation is based. 
Sect. 8 summarizes our findings. 

2. Numerical simulations of quasar microlensing 

The microlensing model used is described in Zackrisson & 
Bergvall ( 2003 1, and is based on the multiplicative magnifica- 
tion approximation (Vietri & Ostriker fl983l l. In this case, the 
magnification, yU to t, due to i microlenses is equal to the product 
of the individual ones: 

/"tot = \\ Hi- (1) 
i 

Although this approximation cannot be expected to repro- 
duce detailed features of microlensing light curves, it has 
been proved useful (e.g. Pei 1993 1 for statistical investigations 
and significantly reduces the calculation time required by ray- 
tracing and similar methods. It has also previously been applied 
(Schneider l 19931 Lacey 1994 > to investigations of quasar light- 
curve statistics. 

When calculating the time-dependent magnification /i(f) re- 
sulting from each compact object, we assume a circular source 
with uniform brightness and the approximate magnification 
formulas for point-mass lenses of Schneider ( 1993 1 and Surpi 
et al. ( 2003 1, relevant for the small and large source regimes 
respectively. The compact objects acting as microlenses are as- 
sumed to be uniformly and randomly distributed throughout the 
simulated observing beam and the influence of external shear 
is neglected. Since the quasars used in this investigation are 
all located in the direction a = 21 h 28 m , 5 = -45°, a univer- 
sal observer velocity will be used. Adopting the velocity of the 
Sun relative to the cosmic microwave background derived by 
Lineweaver et al. (1996 1, the velocity of the observer perpen- 
dicular to the line of sight becomes 307.9 km/s. The velocities 
of quasars perpendicular to the line of sight are assumed to 
be normally distributed with velocity dispersion cr v q S0 = 300 
km/s. 

Unless otherwise stated, all calculations will be carried out 
in the framework of the currently favored A-dominated cos- 
mology, with Om = 0.3, Q.\ = 0.7 and Ho = 65 km/s/Mpc. The 
mean density inside the simulated observing beam is assumed 
to be equal to the mean density of the universe, which implies 
a homogeneity parameter r\ — 1 .0 (Kayser et al. 119971 . 

Zackrisson & Bergvall (2003 1 may be consulted for further 
details on the method used. 

3. Synthetic quasar samples 

In order to be able to draw firm conclusions about possible dis- 
crepancies between simulations and observations, it is impor- 



tant to probe the entire plausible range of microlensing param- 
eter values. When comparing statistics derived from samples of 
observed and simulated light curves, potential selection effects 
must also be considered. 

3.1. The microlensing parameter space 

For each combination of microlensing parameter values listed 
in Tabled 35000 light curves have been generated, uniformly 
distributed in the range of quasar redshifts zqso = 0.13-3.6, 
giving on the order of 15 million light curves. Here, £2 C ompact 
represents the contribution to the cosmological density from 
the compact object population, M compact the mass of the com- 
pact objects responsible for the microlensing effect, cr v compact 
the velocity of compact objects perpendicular to the line of 
sight and Rqso the typical size of the optical-UV continuum- 
emitting region of quasars, normally assumed to correspond to 
some fraction of an accretion disc surrounding a central su- 
permassive black hole. In this investigation, we will assume 
all compact objects to have the same M compact and all quasars 
to have the same Rqso- Even though more complicated mass 
spectra and distributions of source sizes may be present in re- 
ality, we will show that this has no impact on the conclusions 
of this paper. 

We believe that these simulations span the entire microlens- 
ing parameter space relevant for the present situation. The 
cosmological density of compact objects, Q CO mpact, is limited 
upwards by the adopted background cosmology (Am = 0.3, 
Q.A = 0.7), whereas at Q CO mpact < 0.05 the variations induced 
become much too small to explain the observed quasar variabil- 
ity. The time scale of variations induced by lens masses higher 
than M compact = 1 M would typically be longer than the time 
span of present monitoring campaigns, and lens masses lower 
than 10 -5 M Q would give negligible variations except in the 
case of the very smallest Rqso- Where in doubt, single sam- 
ples at 10~ 6 M Q have been generated for test purposes. The 
velocity dispersion of the compact objects is related to their 
preferred whereabouts, i.e. whether mainly located in the field 
or in rich clusters of galaxies. In Zackrisson & Bergvall (2003 1 
we argue, based on results from cold dark matter simulations, 
that cr v compact = 200-600 km/s represents the likely velocity 
span. The parameter Rqso is not well determined observation- 
ally, but in the framework of accretion-disc theory constrained 
to be at least a few times larger than the gravitational radius of 
the central supermassive black hole, which for a 10 8 ~ 9 M ob- 
jects amounts to 1 .5 x 10 n ~ 12 m. The smallest value considered 
here is Rqso = 10 12 m. Values of Rqso significantly larger than 
3 x 10 13 m give very modest light variations and are therefore 
not of interest for the present investigation. 

3.2. Light curves 

Each microlensing light curve generated has a time span of 
25 years, where the time-dependent magnification is estimated 
four times per year within a period of three months. The in- 
trayear variation in magnitudes is then averaged over to give 25 
magnification data points m(t), where m(t) = -2.51og^(f), in 
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Table 1. The set of discrete microlensing parameter values 
which define the grid of synthetic light-curve samples used in 
this study. 



function of a sample of microlensing light curves, representing 
a curve of growth of variability with time lag, is defined as: 
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logfloso ( m ) : 


12 


12.5 


13 


13.5 



close agreement with the procedure used in the observed quasar 
light curves used for comparison. 

Averaging over the intrayear magnification data points has 
the effect of smearing variations that occur on time scales 
smaller than three months and reducing the amplitudes for mi- 
crolensing scenarios where a significant fraction of the varia- 
tions are of this kind. Although the intrayear observations were 
sometimes carried out within a period shorter than assumed 
here, and in some years only one plate was available, tests us- 
ing no magnitude averaging in the yearly data points of the mi- 
crolensing light curves do however show that this discrepancy 
between model and observations has no significant impact on 
the conclusions of this paper. 



In Hawkins ([1996), the variability on time scales shorter 
than one year was studied by taking 16 plates over a period 
of six months. Small short-term variations, supposedly intrin- 
sic in nature, were noted. Even though the present study is fo- 
cused on the long-term variability of quasars, such short-term 
variations will give a small contribution to the variations be- 
tween the yearly averages, and therefore need to be considered 
when assessing the reason for possible discrepancies between 
the observed long-term variability and that predicted by mi- 
crolensing. The observed dispersion among the roughly four 
plates taken each year amounts to 0.13 magnitudes (Hawkins 
2003, private communication), where part of the variability 
comes from the approximately 0.08 magnitude measurement 
error. Assuming the variations due to short-term variability to 
be Gaussian and uncorrected, the effect of both measurement 
error and intrinsic short-term variability can be modelled by 
adding a <r m = 0.065 component of Gaussian noise to each 
yearly m(t). This procedure, which we adopt in the follow- 
ing, is however not entirely consistent for the whole microlens- 
ing parameter space probed, since many combinations of small 
M C o m pact an d small Rqso already predict variability among the 
four intrayear measurements on a level close to that attributed 
to the intrinsic component, and may therefore slightly overpre- 
dict the variability of such scenarios. This will however not 
have any impact on the outcome of the present investigation. 



3.3. Light-curve statistics 

The comparison between observed and simulated light curves 
will be carried out in terms of the average structure function 
and the light curve amplitude. The average first-order structure 



S(J) = 



(2) 



where the sum runs over the N(t) data points for which tj - 1\ = 

T. 

Observationally, the amplitude 6m is defined as the differ- 
ence between the minimum and maximum magnitude observed 
in a quasar when monitored for a certain period of time with 
a particular sampling rate. For the case of variations induced 
purely by microlensing, this reduces to: 



6m = max(m(?)) - min(m(f))- 



(3) 



In the present investigation, the amplitude is derived from 
the first 21 data points of the light curve, in accordance with 
the observational procedure of Hawkins (2000 1. 

3.4. Selection effects 

When comparing the results from microlensing simulations to 
observations, we will use light-curve statistics derived from 
the UVX sample of Hawkins (2000 1, containing 184 objects 
brighter than m B = 21.5 in the redshift range zqso = 0.242 - 
2.21, which has been selected on basis of ultraviolet excess 
(U - B < -0.2). Despite its smaller redshift span, the rele- 
vant selection effects are more easily simulated than in the case 
of variability-selected samples, where the proposed variability 
mechanism must generate approximately the same fraction of 
amplitudes higher than the sample threshold as in real quasars 
in order not to introduce unrealistic biases. As will be shown 
in Sect. 5, microlensing by itself does not seem to be able to 
fulfill this criterion. 

Assuming the colour-criterion of the UVX sample not to 
correlate with variability (see Sect. 7 for a discussion), the most 
important selection effect to consider is the amplification bias. 
A quasar that lacks sufficient intrinsic brightness to be included 
in a flux-limited sample may still be magnified through gravita- 
tional microlensing to reach above the threshold for detection. 
A flux-limited sample will therefore contain an enhanced frac- 
tion of highly magnified objects. As demonstrated in Schneider 
( 1993 ), magnification and amplitude correlate, implying that a 
flux-limited sample should also display an enhanced probabil- 
ity for large amplitudes. This effect may be further boosted by 
the fact that each data point in the observed light curve is also 
associated with an observational error, which statistically in- 
creases the light curve amplitudes. Provided that the error in 
the data point used to define the magnitude of an object is cor- 
related with the amplitude of its light curve, an effect similar 
to that of amplification bias will result from the measurement 
noise alone, increasing the probabilities of high amplitudes af- 
ter the magnitude limit has been imposed. In the present case, 
however, we find that the measurement error is too small and 
the number of data points in the light curve too large for this 
measurement uncertainty bias to have any significant effect on 
the statistical variability properties of the sample. 
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Fig. 1. The observed UVX average structure function (thick 
solid) and the best-fitting (in a least-square sense) average 
structure functions predicted by microlensing scenarios with 
parameter constraints Q CO mpact - 0.15 (thin solid), M compact > 
10 _1 M (thick dashed), and R QS0 > 3 x 10 13 m (thin dashed). 
All simulated average structure functions obeying these param- 
eter constraints but having larger residuals are located below 
the ones plotted. 

The effects of amplification bias is simulated by assigning 
each microlensing light curve a random absolute magnitude 
following the probability predicted by the optical luminosity 
function of Boyle et al. ( 2000 1, converting this into an appar- 
ent magnitude given the background cosmology and assuming 
a power-law continuum (L y oc v~ a ) with slope a = 0.5 for all 
quasars when calculating the ^-correction. The magnification, 
m(0), combined with the corresponding cr m = 0.065 short-term 
variability (measurement error and intrinsic intrayear variabil- 
ity) of the first data point in the light curve is then added to the 
magnitude. Subsequently, all light curves falling on the faint 
side of the UVX magnitude limit are rejected. A selection is 
finally made from the remaining light curves to ensure that the 
redshift and magnitude distributions of the simulated quasar 
sample fully match that of the observed one. 

4. The average structure function 

The observed first order structure function of quasars typically 
consists of three distinct parts (Hughes et al. 119921 . For time 
lags longer than the longest correlation time scale, there is a 
plateau at a value equal to twice the variance of the fluctua- 
tion. For short time lags, there is another plateau at twice the 
variance of the measurement noise. These regions are linked 
by a curve whose slope is related to the nature of the vari- 
ability of the source. In Hawkins ( 2002 ) the observed power- 
law slope (0.20 + 0.01) of the average structure function was 
compared to predictions from three proposed mechanisms of 
variability: The starburst model, accretion-disc instabilities and 
microlensing. The close resemblance to the slope predicted 
by microlensing (0.25 ± 0.03), compared to the much steeper 
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Fig. 2. The observed UVX average structure function (thick 
solid) and the ten best-fitting average structure functions de- 
rived from microlensing scenarios (dashed). 



slopes predicting by the disc instability (0.44 + 0.03) and star- 
burst (0.83 ± 0.08) scenarios appeared to favor microlensing as 
the dominant mechanism. The parameter space of the different 
models was however not probed at any depth. In the case of mi- 
crolensing, structure functions derived from single microlens- 
ing light curves in the literature were used, not average struc- 
ture functions derived from simulated samples. The effects of 
the sample selection bias, the distribution of quasar redshifts 
and the statistical variations between different microlensing 
light curves were thereby neglected. 

Due to the complicated selection procedure of the quasar 
sample from which the observed average structure function 
of Hawkins (2()02i was derived, we will instead compare av- 
erage structure functions derived from our simulated samples 
to an observed average structure function (Hawkins 2003, pri- 
vate communication) derived from the more easily reproduced 
UVX sample. The power-law slope of this average structure 
function is 0.19, very close to the value derived in Hawkins 

Apart from the slope, the scaling of the structure function is 
also important, since this is related to the variance of the light 
curves. It turns out that most of the simulated average structure 
functions lie significantly below the observed one, indicating 
that most parameter combinations generate too little variability. 
If microlensing is to be able to reproduce the observed average 
structure function, certain regions of the microlensing param- 
eter space can therefore be ruled out. In Fig. \l\ we plot the 
observed average structure function in comparison to the av- 
erage structure functions with the highest scaling out of those 
predicted by microlensing in specific parameter ranges. On this 
basis, all scenarios with f2 compact < 0.15, M compact > 10 _1 M G or 
Rqso > 3 x 10°, out of those allowed by the microlensing pa- 
rameter grid defined in Table 1, can be ruled out as inadequate. 
Scenarios with M compact = 10~ 5 M can provide a high degree 
of variability but have a maximum average structure function 
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slope of 0.06 (indicating too much of variability occurring on 
short time scales) and are therefore also unsuitable candidates. 

In Fig. |2 we plot the observed average structure func- 
tion derived from the UVX sample and the ten best-fitting 
(in a least-square sense) simulated average structure functions. 
Even though no single microlensing curve completely repro- 
duces the shape and scaling of the observed average structure 
function, the agreement is reasonable given the limited reso- 
lution of the explored microlensing parameter space. The pa- 
rameters of the best fits lie inside the range Q C ompact = 0.2- 
0.3, M compact = 10 4 -10 3 M , cr v , compact = 200-600 km/s 
and Rqso = 10 12 -10 13 m. This method of determining the 
lens masses required by any viable scenario of microlensing- 
induced long-term quasar variability is in good agreement with 
that of Hawkins ( 1996 1, where the mass of the lensing objects 
was estimated at w 10~ 3 M using simple time-scale argu- 
ments. 

When fitting power-law slopes to the 432 average structure 
functions derived from simulated samples based on the differ- 
ent microlensing parameter configurations of Tabled we find 
the slope to vary in the range 0.004-0.50. The slope of 0.25 de- 
rived in Hawkins ( 2002 ) from single microlensing light curves 
in the literature can therefore not be said to be a robust predic- 
tion of the microlensing mechanism. Instead, the fitted slope is 
quite sensitive to the model parameters assumed. Even in the 
parameter region spanned by the best average structure func- 
tion fits, the range of slopes is 0.06-0.24. 

Even though a large part of the Q CO mpact = 0.2-0.3, 
M compact = 10 4 -10 3 M Q , Rqso = 10 12 -10 13 m microlensing 
parameter space can actually be ruled out on the basis that such 
scenarios significantly overpredict the fraction of high ampli- 
tudes (Zackrisson & Bergvall 2003 1, we make no attempt here 
to determine how close the observed average structure function 
can be reproduced in the absence of these, since this would re- 
quire a significant increase in the parameter resolution. Instead, 
we are content to conclude that there is no strong reason to 
doubt that the observed average structure function can be ex- 
plained by certain combinations of microlensing parameters. 
As we will show, this is however not true for the cumulative 
amplitude distribution or the amplitude-redshift relation. 

5. The cumulative amplitude distribution 

In the absence of any variability-related bias in addition to 
those treated in Sect. 3.4, we should expect microlensing sce- 
narios to approximately reproduce the observed distribution 
of amplitudes in the UVX sample, provided that microlensing 
truly is the dominant mechanism for the long-term optical vari- 
ability of quasars. As shown in Schneider ( 1993 1, the cumula- 
tive amplitude distributions predicted by microlensing simula- 
tions are very sensitive to the model parameters. In Figs.[3]and 
|4]we show examples of the cumulative amplitude distributions 
predicted for different lens masses in the case of Rqso = 10 12 m 
and 3 x 10 13 m, respectively. 

In Hawkins ( 1996 1, it was argued that the amplitude distri- 
butions predicted by the microlensing simulations of Schneider 
(1993 1, which assumed Rqso = 10 13 m and an Einstein-de 
Sitter Universe, could probably provide a good match with the 
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Fig. 3. The cumulative probabilities P(< 5m) of observing 
quasars with amplitudes lower than 5m for the UVX sam- 
ple (thick solid). The other lines represent the corresponding 
distributions derived from microlensing simulations assuming 
^compact = 0.3, Rqso = 10 12 m, cr v compact = 400 km/s and 
^compact = 1 M e (thin solid), 10 1 M (thick dashed), 10~ 2 M 
(thin dashed), 10 3 M (thick dash-dotted), 10 4 M Q (thin 
dash-dotted), 10~ 5 M Q (thick dotted). 

observed distribution of amplitudes for M compact = 10~ 3 M , 
provided that Rqso = 3 x 10° m instead. Our model confirms 
this to be the case for the Einstein-de Sitter cosmology, but not, 
as we show in Fig.0] for the Dm = 0.3, Qa = 0.7 cosmology 
assumed here. At such large Rqso, no microlensing scenario 
even comes close to reproducing the observed amplitude dis- 
tribution. This is in good agreement with our conclusions from 
the analysis of average structure functions, where scenarios as- 
suming Rqso ^ 3 x 10 I3 m were also ruled out on the basis of 
giving to little variability. 

In Fig. |3] many of the cumulative probability distributions 
predicted by microlensing scenarios can be seen to clearly cross 
over to the right-hand side of the corresponding observed dis- 
tribution, indicating an overprediction of the fraction of the 
very highest amplitude objects. This fact was used in Schneider 
( 1993 1 and Zackrisson & Bergvall ( 2003 1 to infer upper limits 
on ^compact for different M compact . Here, we will instead concern 
ourselves with another discrepancy between the observed and 
predicted amplitude distributions more directly related to the 
question of whether microlensing really is responsible for the 
observed light variations or not. 

The selection criterion of 5m > 0.35 has proved useful 
(Hawkins 2000 1 for building observational samples of quasars 
extending to redshifts higher than the zqso ~ 2.2 limit of sam- 
ples based only on ultraviolet excess. The fraction of objects 
with amplitudes higher than 0.35 magnitudes predicted by the 
dominating mechanism of variability therefore has important 
implications, since it is related to the completeness of such 
variability-selected samples. If an adequate fraction of ampli- 
tudes higher than 0.35 magnitudes cannot easily be reproduced 
by microlensing, there is no way to simulate the properties 
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Fig. 4. Same as Fig.|3]except Rq 5Q = 3 x 10 13 m. 



of variability-selected samples using the 6m > 0.35 criterion 
(like in the VAR sample of Hawkins 2000 ), since this would 
mean rejecting an unrealistically large fraction of simulated 
light curves. 

As can be seen in Figs.|3]and |4] none of the microlensing 
scenarios explored manage to reproduce the observed fraction 
of amplitudes higher than 0.35 magnitudes. The microlensing 
simulations instead predict large numbers of objects below this 
limit. This turns out to be a general feature of our simulations, 
and not just of the small number of parameter combinations 
used in these plots. It should be noted that the absence of am- 
plitudes smaller than 0.2 magnitudes in the predicted distribu- 
tions is almost entirely due to the presence of simulated obser- 
vational errors and intrinsic intrayear variability in each data 
point of the light curves. Without the added short-term noise, 
significant fractions of even lower amplitudes would be pre- 
dicted. 

In Fig.[5]we compare the observed fraction of objects with 
amplitudes higher than 0.35 magnitudes in the UVX sample 
with the corresponding microlensing predictions as a func- 
tion of M compact and R QS0 in the case of Q C om P act = 0.3 and 
o"v,compact = 600 km/s (the scenario giving the most variability). 
In no case can microlensing fully reproduce the observed frac- 
tion. Given the trend of variability with lens mass in the case 
°f ^qso = 10 12 m, one could however imagine that extending 
the investigation to lower lens masses than 10~ 5 M Q might im- 
prove the agreement, but tests show that the predicted fraction 
of amplitudes above 0.35 magnitudes does not become signifi- 
cantly higher at 10~ 6 M Q . Such scenarios can also be ruled out 
on the basis of giving too much short-term variability (structure 
function slopes much flatter than observed). 

Hawkins ( 2002 1 shows that the Seyfert galaxies in his sam- 
ple are dominated by a different kind of variability. Since the 
UVX sample used here also contains a number of such ob- 
jects, one might expect that this could affect the comparison. 
However, following the procedure of Hawkins (2002i of re- 
moving objects with zqso < 0.5 and Mb > -23 to get a clean 
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Fig. 5. The fraction of objects with amplitudes higher than 0.35 
magnitudes predicted by microlensing scenarios for different 
M compact , assuming Q comp act = 0.3, cr ViCompact = 600 km/s and 
R QS0 = 10 12 m (thick solid), 3 x 10 12 m (thin solid), 10° m 
(thick dash-dotted), 3 x 10 13 m (thin dash-dotted). The thick 
dashed line marks the observed fraction (93%) of the UVX 
quasar sample. 

sample of quasars does not improve the agreement, but actually 
makes it slightly worse. 

We conclude that the fraction of amplitudes higher than 
0.35 magnitudes cannot fully be reproduced by our present 
microlensing simulations. The gap between observations and 
predictions is however quite small for lenses in the range 
M compact = 10~ 5 -10~ 4 M , as seen in Fig.0 

6. The amplitude-redshift relation 

In Hawkins l|2000l, the observational relation between mean 
amplitude and redshift was explored. In the UVX sample, the 
mean amplitude was found to increase up to a redshift around 
one, followed by a roughly constant mean amplitude towards 
higher redshift. This was attributed to a magnitude effect in the 
sense of higher luminosity quasar having smaller amplitudes. 

The most naive expectation from microlensing scenarios 
would perhaps be a mean amplitude which continuously in- 
creases with redshift due to increased optical depth to lens- 
ing. Qualitatively, the observed behaviour of the mean ampli- 
tude may however be explained by microlensing-induced vari- 
ations, provided that the size of the continuum-emitting region 
at the rest wavelength observed is larger for the more luminous 
quasars found at high redshifts. 

In Fig. [6] the empirical amplitude-redshift relation is plot- 
ted together with the corresponding relations derived from the 
432 different microlensing parameter configurations allowed 
by Table □ Since every microlensing curve assumes Rqso to 
be constant with redshift, the fact that the shapes of the simu- 
lated and observed relations are very different does not rule out 
microlensing as a plausible mechanism. More serious is how- 
ever the fact that no microlensing scenario appears capable of 
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Fig. 6. The empirical relation between mean amplitude and red- 
shift for the UVX quasar sample (thick black) and the corre- 
sponding relations derived from the 432 different microlensing 
parameter combinations generated (thin grey). The error bars 
indicate the uncertainty in the position of the observed mean. 

explaining the high mean amplitudes observed at zqso < 1, 
where one third of the quasars in the UVX sample are located. 
Even though relaxing the approximation of constant Rqso may, 
in principle, create a plateau in the amplitude-redshift relation, 
the mean amplitudes can never lie higher than the maximum 
microlensing prediction of Fig. |5J provided that microlensing 
really is the dominating mechanism for variability. This dis- 
crepancy between observations and simulations is much larger 
than that discovered in connection to the cumulative amplitude 
distributions. We therefore conclude that, at least for quasars at 
Zqso < 1, significant contributions to the long-term variability 
must come from a mechanism other than microlensing in order 
to explain the observed amplitude-redshift relation. 

7. Discussion 

Even though our microlensing simulations appear capable of 
reproducing the observed average structure function of quasars, 
they fail to explain both the observed fraction of objects with 
amplitudes higher than 0.35 magnitudes and the high mean am- 
plitudes observed at zqso < 1 ■ This does not rule out the pos- 
sibility that microlensing may still contribute to the variability 
at some level, but indicates that another mechanism must also 
contribute significantly. 

Since any mechanism other than microlensing contributing 
to the observed variability would have to give largely achro- 
matic and symmetrical variations to be consistent with obser- 
vations, this means that microlensing cannot be the only mech- 
anism involved possessing these properties. The exciting possi- 
bility of deriving properties of dark matter by comparing results 
from microlensing simulations to the observed properties of 
quasar light curves is thereby severely obstructed, since any mi- 
crolensing signal that may be present in the light-curve statis- 
tics is likely to be blended with a different kind of variability 



with similar characteristics. Although the chromaticity of pro- 
posed mechanisms of intrinsic variability like accretion-disc in- 
stabilities and supernova explosions have not been thoroughly 
investigated, both of these mechanisms appear capable of gen- 
erating statistically symmetric variations on long time scales 
(e.g. Kawaguchi et al. 1998 Aretxaga 1997 1. The one piece 
of evidence pointing uniquely in the direction of an extrinsic 
mechanism is however the lack of cosmological time dilation in 
the characteristic time scales of the variations (Hawkins 2001 1. 
Provided that no efficient extrinsic variability mechanism other 
than microlensing can be found, one possible way of setting a 
lower limit to the contribution from microlensing could there- 
fore be to investigate how large the microlensing population 
would have to be to eradicate the cosmological time dilation 
signal believed to be a feature of the intrinsic variability com- 
ponent. 

How robust are then the results presented here? So far 
we have considered only populations of compact objects and 
quasars characterized by single values of Rqso and M compact . In 
reality, extended distributions of quasar radii (possibly corre- 
lated with quasar luminosities) and lens masses may of course 
be present. In the framework of the microlensing model used, 
such distributions pose no threat to our conclusions, as long as 
our simulations truly probe the entire range of plausible param- 
eter values. Since no single parameter combination explored 
can reproduce the high mean amplitudes observed at low red- 
shift, neither will any mixture. 

The results presented assume Qm = 0.3, Qa = 0.7, but the 
conclusions are not sensitive to minor variations in the adopted 
cosmology. In Fig.0 we present the predicted mean amplitude- 
redshift relation in the case of Q CO mpact = Qm, M compact 
1 : M B , Rqso = 10 12 m, cr vcompact = 400 km/s (one of the 
scenarios giving the highest mean amplitude in the lowest red- 
shift bin of Fig. [6jl for a number of different combinations of 
Qm and Qa- Clearly, variations like (Qm = 0.2, Qa = 0.8) 
and (Qm = 0.4, Qa = 0.6) allowed by current uncertainties 
in these parameters do not provide an obvious way around the 
discrepancies between simulations and observations. Only in a 
cosmology with a matter density significantly higher, like in the 
case of an Einstein-de Sitter Universe (Qm = 1.0, Qa = 0.0), 
is there a reasonable agreement at the lowest redshifts. 

Our simulations also assume rj = 1, which is admittedly 
not self-consistent for the entire microlensing parameter space 
investigated. This parameter describes the fraction of homoge- 
neously distributed matter, where the scale of clumps in the re- 
maining matter per definition is on the same order of magnitude 
as the angles involved (Kayser et al. 1997 ). The rj parameter af- 
fects the angular diameter distances going into the magnifica- 
tion contribution from each lens and the luminosity distances 
used to make the conversion from absolute to apparent magni- 
tudes. When the microlensing source is large and the variations 
are due to Poisson fluctuations in the number of compact ob- 
jects in front of the source (Surpi et al. 2003 1, rather than due to 
rare and isolated microlensing events, the mean density of the 
simulated observing beam becomes approximately equal to the 
mean density of the universe, which makes rj = 1 reasonable. 
However, when the angular extent of the source is small com- 
pared to that of the lenses, the density in the simulated observ- 
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Fig. 7. The empirical relation between mean amplitude and red- 
shift for the UVX quasar sample (thick solid) and the corre- 
sponding predictions from microlensing simulations assuming 
different background cosmologies: The standard Dm = 0.3, 
Oa = 0.7 case (thin solid), Qm = 0.4, Qa = 0.6 (thick dashed), 
Q M = 0.2, Q A = 0.8 (thin dashed) and Q M = 1.0, Q A = 0.0 
(thick dash-dotted). All simulations assume Q C ompact = Qm, 
M compact = 10~ 4 M Q , R QS0 = 10 12 m and o- v , compact = 400 
km/s. The error bars indicate the uncertainty in the position of 
the observed mean. 

ing beam varies violently with time. In this case tj < 1 is more 
appropriate. This situation arises naturally when modelling the 
gravitational lensing of supernovae (e.g. Mortsell 2002 1, which 
are generally assumed to be less extended than quasars, but 
also in the present investigation when small values of Rqso 
are considered. By testing different values of rj in the range 
j] = 0-1, we have however found that variations in this param- 
eter have only a negligible impact on the amplitude distribu- 
tion and the mean amplitude-redshift relations predicted by our 
simulations, and therefore pose no threat to our conclusions. 

A more serious concern than the assumed values of our 
model parameters is a number of simplifying approximations 
on which the microlensing simulations are based. 

The multiplicative magnification approximation, which de- 
spite having been used in similar contexts several times be- 
fore (e.g. Schneider ll993l Lacey[l994i and having been shown 
to give reasonable accuracy when it comes to the statistics of 
magnification distributions (Pei 1993 1, has not been thoroughly 
tested for the case of light-curve statistics. A particular con- 
cern is that it may underpredict the amplitudes, since the high- 
amplification features associated with caustic crossings are not 
well-reproduced in this approximation. This can e.g. be seen by 
comparing the light curves in Schneider ( 1993 1 where the mul- 
tiplicative magnification approximation is used, with those in 
Lewis et al. ([1993 1, which rely on a ray-shooting approach. The 
naive expectation would be for inclusion of light-curve features 
associated with caustic crossings to increase the mean ampli- 
tudes and flatten the average structure function due to introduc- 
tion of more short-term variability. The actual effect will how- 



ever depend on the details of the sampling of the light curves 
and is hard to quantify without making a direct comparison be- 
tween the two methods. If a proper treatment of caustic cross- 
ings would in fact give a significant contribution to the pre- 
dicted amplitudes, this would also strengthen the upper limits 
on O compact inferred by Zackrisson & Bergvall ( 2003 1, where a 
large part of the relevant small source microlensing parameter 
space has already been ruled out. For this reason, it is doubtful 
whether relaxing the approximation of multiplicative magnifi- 
cation really can bring the model into better agreement with 
observations. 

Our simulations furthermore assume all compact objects to 
be randomly and uniformly distributed with constant comoving 
density throughout the simulated observing beam, an approach 
commonly referred to as the Press & Gunn dl973> approxima- 
tion. Wyithe & Turner (2002 1 showed that although this ap- 
proximation can give very misleading results if only known 
stellar populations contribute to the microlensing effect, it is 
useful for estimating the microlensing optical depth in the case 
where a significant fraction of the dark matter is in the form of 
compact objects, as in the scenario explored here. This does not 
necessarily imply that this approximation is valid when gener- 
ating light-curve statistics, since correlated lens positions could 
give rise to correlated magnifications in the simulated time se- 
ries and therefore affect the distribution of amplitudes. When 
evaluating such possible effects, one should instead of the iso- 
lated dark matter halos assumed in Wyithe & Turner more re- 
alistically use the mass distributions predicted when sending 
the observing beam through cosmological N-body simulations 
of large-scale structure. The naive expectation for the inclusion 
of clustered lenses is however an increased bimodality in the 
distribution of amplitudes, as some simulated observing beams 
will contain huge concentrations of compact objects whereas 
other will be essentially devoid of matter. This is likely to de- 
crease the fraction of light curves with amplitudes above 0.35 
magnitudes, since a growing number of quasars should display 
practically no microlensing variability. The demonstrated dis- 
crepancy between microlensing predictions and observations in 
Fig.|5]are therefore expected to grow even further as effects of 
clustering are introduced. 

Testing the reliability of these approximations is not only 
important for investigations of the present kind, but also for 
modelling the time-dependent microlensing magnification of 
other light sources at cosmological distances, like supernovae 
and the afterglows of gamma-ray bursts. To conduct the re- 
quired investigation does however represent a huge undertaking 
and is not within the scope of this article. Instead, detailed tests 
of the multiplicative magnification approximation and the Press 
& Gunn approximation in the context of light-curve statistics 
are intended to form the basis of a future paper. 

Assuming that our simulations really are sufficiently reli- 
able, could it then be that the discrepancy between microlens- 
ing predictions and observations is simply due to a miss- 
ing fraction of low amplitude quasars in the UVX sample? 
Ultraviolet excess is normally not assumed to correlate with 
variability (e.g. Hawkins & Veron fT995l . but provided that the 
quasar continuum radiation originates from a large accretion 
disc which becomes bluer towards the center, one could still 
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argue that a lensing event with small impact parameter would 
statistically both increase the amplitude and decrease U - B, 
thereby introducing a correlation between these two parame- 
ters. There are however two reasons why this effect is unlikely 
to produce a bias of the required strength. Firstly, because the 
observed achromaticity of the quasar light curves are inconsis- 
tent with strong colour gradients across the accretion disc if the 
variations truly are caused by microlensing in the large source 
regime, and secondly because huge numbers of quasars would 
have to be present on the redward side of the U - B < -0.2 
threshold, which does not appear to be the case (e.g. Veron 

8. Summary 

Using a microlensing model based on the multiplicative mag- 
nification approximation and the assumption of uniformly and 
randomly distributed compact objects, we find that microlens- 
ing can reasonably well reproduce the observed average struc- 
ture function of quasars, but fails to reproduce the high fraction 
of objects with amplitudes higher than 0.35 magnitudes and the 
mean amplitudes observed at zqso < 1 ■ This indicates that an- 
other mechanism must also contribute significantly to the long- 
term optical variability and severely complicates the task of us- 
ing light-curve statistics from quasars which are not multiply 
imaged to isolate properties of any cosmologically significant 
population of compact objects which may be present. The one 
piece of evidence pointing uniquely in the direction of an ex- 
trinsic variability mechanism is however the lack of cosmolog- 
ical time dilation in the characteristic time scales of the varia- 
tions. If no other efficient mechanism for extrinsic variability 
than microlensing can be found, one possible way of setting a 
lower limit to the contribution from microlensing could there- 
fore be to investigate how large the microlensing population 
would have to be to eradicate the cosmological time dilation 
signal believed to be a feature of any intrinsic component of 
variability. 
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